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Abstract 



Noncolliding diffusion processes reported in the present paper are A'^- 
particle systems of diffusion processes in one-dimension, which are conditioned 
so that all particles start from the origin and never collide with each other 
in a finite time interval (0, T), < T < oo. We consider four temporally 
inhomogeneous processes with duration T, the noncolliding Brownian bridge, 
the noncolliding Brownian motion, the noncolliding three-dimensional Bessel 
' bridge, and the noncolliding Brownian meander. Their particle distributions 

at each time t G [0, T] are related to the eigenvalue distributions of random 



, matrices in Gaussian ensembles and in some two-matrix models. Extreme 

\ values of paths in [0, T] are studied for these noncolliding diffusion processes 

and determinantal and pfaffian representations are given for the distribution 



functions. The entries of the determinants and pfaffians are expressed using 
Q \ special functions. 

KEY WORDS: Noncolliding diffusion processes, random matrix theory, ex- 
treme value distributions, determinants and pfaffians, special functions ^ 
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5^ ; 1 Introduction 

For N >2, consider the following region in M^, 

= ja; = (xi,X2, . . .,Xn) e : Xi < X2 < ■ ■ ■ < Xn^ (l-^) 



which is called the Weyl chamber of type A^r-i in the representation theory [9]. 
By the Karlin-McGregor formula [10], the density at y E of an A^-dimensional 
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Brownian motion at time t > 0, which starts from x e at time 0, and is 
restricted on the event that it stays in W;^ during a time interval [0,t\, is given by 

f^{t,y\x) = ^ det p{t,yi\xj) 



where p{t^y\x) is the heat kernel, p{t,y\x) = e /{2t) j^^ixt. It can be regarded 

as the transition probability density of the absorbing Brownian motion in from 
a; e to t/ e with duration t > 0. 

Set < T < oo, a = (ai, . . . , a^r), b = . . . , ftjv) e W^. Then we consider the 
A^-particle system of one-dimensional Brownian motions starting from a configura- 
tion a at time and arriving at the configuration b at time T, conditioned never 
to collide with each other in [0,T]. The probability density at the configuration 
X e at time t e [0, T] is given by 







a) 




a) 



PN,T{t^ ^; a, b) = M^^KiV (1-2) 



We will call the above mentioned process the noncolliding Brownian motion from a 
to h with duration T and write it as X'^'^it) = (Xf'^(t), .... X^^{t)),t E [0,T]. 
We can show that the limit of (1.2) in a ^ = (0, 0, . . . , 0), b ^ 0, is given by 

, lim p^,t(^, a, b) = a;; 0, 0) (1.3) 

|a|,|bKO 

= g^uE^a^jMt)^), te[0,T], o^eW^ 



with ar(t) = v^t(l - t/T) (Proposition 13 in [15]). Here q'^^^{x;a^) denotes the 
unitary ensemble (GUE) with variance cx^, 



probability density of eigenvalues x G of random matrices in the Gaussian 



where |a;p = '^f^ix'j and r(A:) is the Gamma function. Note that, when k E 
N = {1, 2, . . . }, r(/c) = (A; — 1)!. We regard (1.3) as the probability density at a 
configuration x E at time t E [0, T] of the noncolliding Brownian bridges with 
duration T denoted by X°'"(i) = {X^'°{t), X^^{t)),t E [0, T]. 

When we make the configuration b at time T be arbitrary in W^, we have 
another temporally inhomogeneous system of noncolliding Brownian motion, which 
is denoted by X^'^(t) = (Xf'*(t), . . .,X^^{t)),t E [0,T]. The probability density 
at a; e at time t E [0, T] is given by 

(f ^.:,J,^- ^NiT-t,x)f^{t,x\^) 

p^T[t,x,a.,K) ./-A/^ N ' *e[U, JJ, 

J^nK-'- > a; 
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where 



J\f^{s,x)= / dyf^{s,y\x), s > 0, a; e W^, 



N 



is the probabihty that the absorbing Brownian motion in W;^ starting from £C G W 
is not yet absorbed at any boundary of the region and is surviving inside of it 
at time s > 0. For an even integer n and an antisymmetric n x n matrix A = {aij) 
we put 



Pf 

l<ij<n 



(»/2)! 



E 



a:a(2k-l)<a(2k),l<k<n/2 



where the summation is extended over all permutations a of {1,2, ... ,n) with re- 
striction a{2k — 1) < <y[2k), k = 1, 2, .... n/2. This expression is known as a pfaffian 
(see, for example, [24]). By using the de Bruijn identity [5], which will be given as 
Lemma 3.1 in Section 3, we have the formula [12, 13] 



X^(s,x)^ Pf 



l<i,j<N 



X 



where 



N ^ 



N, 

N + 1, 



x/y/2s — {xi/ \/2s, . . . , xn/\/2s), and 




with 



a N = even, 
if = odd. 



if l<^,J<iV, 

if 1 < i < iV, j = iV + 1, 

iii^ N + l,l<j < N, 



e-" dv. 



:i.4) 



In [12, 13], the process = lim|a|^o ^"'^W = (Xf '^(t), . . . , X°'^(t)), t e 

[0,T] is constructed, where the probability density at a; e at time t e [0,^] is 
given by 



rpN{N-l)/4^-N^/2 



II {xj-Xi)Af{T-t,x). (1.5) 



l<i<j<N 



It has been shown that (1.5) exhibits a transition from the eigenvalue distribution of 
GUE to the eigenvalue distribution of another ensemble of random matrices called 
the Gaussian orthogonal ensemble (GOE) as t — )■ T [13, 14]. In other words, there 
establishes an interesting correspondence [12] between the temporally inhomoge- 
neous system of noncoUiding Brownian motion X^'^(t),t G [0, T] and the two-matrix 
model of Pandey and Mehta [21] in the random matrix theory [20]. 
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The equivalence of the distribution of and the eigenvalue distribution of 

random matrices in GUE with variance aT{tY for t e [0,7] shown by Eq.(1.3), and 
the correspondence between X^'^{t),t G [0, T] and the two-matrix model of Pandey 
and Mehta are very useful to perform computer simulations of the conditional diffu- 
sion processes and e [0,T] [19]. Figure 1 shows samples of paths 
generated by computer simulations for these two processes. 




Time Time 



(a) (b) 



Figure 1: Samples of paths for (a) and (b) X°''^(t),t G [0,T], generated by 

simulating the corresponding eigenvalue processes of random-matrix models. 

In the present paper, we study the distributions of the extreme values defined 

by 

LHN,T) = mm XHt) = mm xUt), (1.6) 

R^(N,T) = max X^(t) = max xUt), 

te[o,r] te[o,T] 

for H = (a, b) G X W^, = (0, 0), and (0, M). 

In this paper, we also consider the noncolliding processes "with a wall at the 
origin". Instead of the Weyl chamber of type A^v-i given by (1.1), we consider the 
Weyl chamber of type Cjv, 



W 



|a; — {xi, X2, . . . , xn) G : < Xi < X2 < ■ ■ ■ < xn^. 



The density &t y E of an A/'-dimensional Brownian motion at time t > 0, which 
starts from x G at time 0, and is restricted on the event that it stays in 
during a time interval [0, t], is given by the Karlin-McGregor formula as 

fS{t,y\x) = Jet \Af{t,y^\x,)], (1.7) 

l<t,j<N L 

where P^^]^^\t,y\x) is the transition probability density of the one- dimensional ab- 
sorbing Brownian motion with an absorbing wall at the origin. By the reflection 
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principle of Brownian motion, it is given as 

p'fhr\'t, y\x) = p{t, y\x) - p{t, y\ - x) 

^ ^ f -(x-y)V{2t) _ g-(a;+j/)V(2t)]. 

Eq.(1.7) gives the transition probability density of the absorbing Brownian motion 
in from x to y with duration t > 0. 

For 0<T<oo, a, be W^, we consider the A^-particle system of one- 
dimensional Brownian motions starting from a configuration a at time and arriving 
at a configuration b at time T, which is conditioned so that there is no collision be- 
tween any pair of particles and that all particles stay positive in a time interval 
[0, T]. The probabihty density at a; e at time t e [0,T] is given by 

The Brownian motion conditioned to stay positive for a finite time interval [0, T] is 
called the Brownian meander with duration T [22] . So we call the A^-particle system, 
whose probability density is given by (1.8), the noncoUiding Brownian meander from 
a to b with duration T [16] and write it as V'^it) = {Y-^'^{t), Y^'^{t)),t e [0, T]. 

The three-dimensional Bcssel bridge with duration T is the conditional Brow- 
nian motion such that it starts from the origin, stays positive in t e (0,T), and 
first returns to the origin at time T. Then the A'"-particle system obtained by the 
limit Y^'°{t) = lim|a|,|bKo^^^'''(0 = {Y°'°{t), . . . ,Y°'^{t)),t e [0,r] can be called 
the noncoUiding three-dimensional Bessel bridge. It was called the system of nonin- 
tersecting Brownian excursions in [25]. The probability density at a? e at time 
i e [0, T] is given by 

p'^^Jt,x-Q,Q) = q%{x-aT{tf) 



with (Jr(t) = \/t{l — t/T) (Proposition 14 in [15]). Here q%{x] cr^) denotes the prob- 
ability density of positive eigenvalues x e of random matrices in the ensemble 
called the class C with variance cr^ [1], 



-7V(2Ar+l) ^ 



When we make the configuration b at time T be arbitrary in W^, we have 
another noncoUiding Brownian meander, which is denoted by 

Y^'^+{t) = {Y^'^+{t),...,Y^^+{t)), t e [0,T] with M+ = {x e M : X > 0}. The 
probability density at cc e at time t e [0,T] is given by 

a f.^.^E) ^fSiT-t,x)f^{t,x\^) 
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where 



M^{s,x)= [ dyf^{s,y\x), s > 0, xeW%, 



is the probabihty that the absorbing Brownian motion in starting from x E W% 
is not yet absorbed at any boundary of and is surviving inside of it at time s > 0. 
We can prove the formula [12, 13] 



J\f^(s, x) = Pf 



l<i,j<N 



X 



where x/ v2s = {xi/y/2s, . . . , xn/V^s) 



'^{yi,yj), 

0, 



if l<^,J<iV, 
ifl <i < NJ = N +1, 
iii = N + l,l<j <N, 
if i = j = AT + 1, 



with (1.4) and 



^'(«l,M2) 



2 

TT 



ni PU2—UI 

J U\—U2 

U2 rui+U2 

dvi / d^;2e-"?-("i-"^)' 

Ml J U2—U1 



The process = hmia^o 1""'''+ (t) = (Y.^'^^t), . . . ,Y^^^^ {t)),t E [0,T] is 

studied in [13, 18], where the probabihty density at a; e at time t e [0,T] is 
given by 



rpN"^ /2^-N{2N+l)/2 



:i.9) 



N 



l<i<j<N 



k=l 



It has been shown that (1.9) exhibits a transition from the eigenvalue distribution 
of random matrices in the ensemble called the class C to the eigenvalue distribution 
of the class CI as t ^ T [18, 15]. 

We can perform computer simulations of a matrix-valued Brownian bridge such 
that its distribution is the same as the distribution of random matrices in the class 
C with variance axit) = \/ t{l — t/T), t G [0,T], and a matrix- valued Brownian 
motion, whose distribution at time t e [0, T] changes continuously from the class 
C distribution to the class CI distribution of random matrices as t ^ T. By nu- 
merically calculating eigenvalues of these two matrix-valued diffusion processes in 
t e [0, T], we can draw samples of paths of noncolliding diffusion particles for Y^'^(t) 
and y"'^+(t),t e [0,r], as shown in Figure 2 [19]. 
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(a) 



(b) 



Figure 2: Samples of paths for (a) and (b) Y^'^+{t),te [0,T]. 

For Y^, = (a, b) e x W^, (0, 0), and (0, R+), we will report the distribu- 
tions of the maximum value 

hHN,T) = max YHt) = max vlU). (1.10) 

te[o,T] te[o,T] 

2 Results 

For /c e No = {0, 1, 2, . . . }, let Hk{x) be the k-ih Hermite polynomial 

-l)™(2x) 



m=0 



■m\{k — 2m)\ ' 



where [a] denotes the largest integer that is not greater than a. We define the 
function of ti, v G M with an index k eNq, 

ek(u,v) = ^//fe(iin + ^;)e-(""+'')'. (2.1) 
If we consider the following version of Jacobi's theta function 

Imy > 0, where we have set z — eJ^^J~^^ and q — e^^^'^ ^ the reciprocal relation 

1/2 



holds (sec Sec. 10. 12 in [2], Sec.A.3.1 in [8]). Combining this functional equation 
with the fact 



OO h. 
2 ... s M 



fe=0 

7 



the following identities are derived for e No, r/, ^ G M (Lemma 3 in [19]); 



exp 



TT 



-n 



The following two theorems are our main results. 
Theorem 2.1 For i,r > 0, 

p(^-e < L^'^{N,T),R^'^{N,T) <r) 



-1) 



N 



2N{N-l)/2 Yl^^^ T{k) 



X det 

l<i,j<N 



e 



i+j-2 



2(£ + r) 2^ 



2T V2T 



V V2r 



and 



P{ - £< L°'^(Ar, T), i?°''^(Ar, T) < r) 

(_l)iV(^^+l)/2 



2iv(iv-i)/4]-f^^^r(A;/2) i<i,j<N 



Pf 



where 



G 




if l<i<N,j = N + l, 
if i = N + l,l<j < N, 
0, if i = j^N + l, 



with 



2{i + r) 2£ 



2T V2T 



+ a; +(-iye,_i 



2(£ + r) 



-,a; 



7^ = 



e/V2T 



zf^(x)dx, 



det 



-£/V2T<xi<a;2<r/V2T 

Theorem 2.2 For /i > 0, 

(-1)^ 



2;f(,xi) 2;f(.X2) 



dxidx2- 



det 



2^' rifeLi r(2A;) i<i.^<^ 



e 



2(i+i-l) 



2/1 



,0 



and 



Pf 



2^M/2nl;=ir(A;) i<^,,<iv 
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where 



G 



h 



if 1 <i,J < N, 

I?, if l<i<N.] = N + l, 

-If, if ^ = iV+l,l<J<A^, 

0, if i = j = N+l, 



with 



2h 



2T 



X 



zl'{x)dx, 



rP= [ 

Jo- 



det 



0<xi<X2<h/V2T 



zf{xi) zf{x2) 



dxidx2- 



Remarks. 



(1) By the scaling property of Brownian motion, we have the equahties in dis- 
tribution 

L^{N, T) = VtL^{N, 1), R^{N, T) = VfR^{N, 1), 
for tl = (0,0),(0,M), and 

h\n,t) = Vth\n,i) 

for jl = (0,0), (0,R+), VT > 0. As a matter of fact, the probabihty distribu- 
tions (2.3) and (2.4) in Theorem 2.1 are functions of 1/\/T and r/\/T, and 
(2.5) and (2.6) are of h/\/T, respectively. 

(2) For = 1, (2.5) gives 

p(ff°.»(l,T)</.) = -IY^hJ^u) 



Ah'^n^ 



since H2{x) — Ax^ — 2. This is a classical result for the height distribution of 
the three-dimensional Bessel bridge. It should be remarked that this result for 
distribution is equivalent with the fact on moments 



E 



H' 



°-°(l,T))' 



TT 



rp\ m/2 



with 



e(m) = im(m - l)7r-/^r (y) C(m), 



where C(^) is the Riemann zeta function 

^H = J];^, Rem>l. 

neN 

See [3] for interesting relations of the probability laws of one-dimensional con- 
ditional Brownian motions to the Jacobi theta functions and the Riemann zeta 
function. 

(3) In [11], N = 2 case of the formula (2.5) of Theorem 2.2 and moments cal- 
culated from it have been intensively studied, which are related to the double 
Dirichlet series. 

(4) The formula (2.3) of Theorem 2.1 gives the joint distribution of two extreme 
values L^'^{N,T) and i?0'0(Ar,T) of the process X°'"(i),i e [0,T]. The distri- 
butions of single variable R^'°{N,T) was studied in [7] and [23]. Feierl [7] also 
reported the distribution of the width (or "range" ) defined by 

W^'°{N, T) = R^^^{N, T) - L^'^{N, T) 

= max (x^/{t)-X^,^\t)). 

(5) A part of results recently reported by Borodin et al. [4] may be stated as 
follows in the present notations; if is even, for each fixed T > 0, 

max , ^ X^'°(t)4y"'°(r/2), 
*6[0,T/2] ^2(1 - t/T) ^ ^ N/2\ I )^ 

where Y^'^^{T/2) in the RHS denotes the position of the top particle at time 
T/2 of the iV/2-particle system l"°'°(t) = {Y^'''{t), . . . ,Y^f^{t)),t e [0,T]. 
The distribution of the top path Y^'l^^{t),t e [0,T] was studied by Tracy and 
Widom [25]. 

(6) The distribution of the "maximum height" H^'^{N, T) of the process Y^'^{i), t G 
[0, T] has been studied by Feierl [6] and by Schehr et al. [23] in different context 
and by different methods. See [19] for comparison. 



3 Proof of Theorem 2.1 

Let be the indicator function of a condition ui. Assume that — ^ < ai, — ^ < 

61, ajv <r,hN < r. By definition of the process and (1.6), for (a, b) e W;^ x W;^ 

P(^-£ < L^'^{N,T),R'''^{N,T) < r) (3.1) 
= e(i{-£ < X^'''{t),X^/{t) < r, e [o,r]}) 
= < X^'^it) <r,l<yi<N,yte [0,T]}). 
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For £,r > 0, we introduce the following restricted region in W^, 

W^{-i, r) = {a; = {xi, . . . , x^) e : < Xi, x^v < r}. 



Then we consider the absorbing Brownian motion in this region starting from a G 
W^(— £, r). Let the transition probability density of a one-dimensional Brownian 
motion in an interval {—£, r), where two absorbing walls are put at x — —£ and x — r, 

^® ^'abs' V'l' r)>^ ^ 0- Then, by the Karlin-McGregor formula, the probability 
density that the absorbing Brownian motion in W^(— £, r) arrives at b G W^(— £, r) 
at time T > 0, avoided from any absorption at boundary of the region W;^(— £, r), 
is given by 



/^(-^-nT,b|a)= det p™(T,6,|a,) 

l<i,j<N L 

The probability density p^^^'^\t,y\x) is obtained by solving the diffusion equation 
du{t,y)/dt = {l/2)d'^u{t,y)/dy'^, —i<y<r with the Dirichlet boundary condition 
u{- , —i) = u{- , r) = and with the initial condition u{0, y) — Sx{y), —£< x,y < r. 
The solution is given by 

t > 0,x,y e {-i,r). 

Then the following equality will be estabhshed, 

e(i{-£ < Xf'^t) < r, 1 < Vi < TV, G [0,T]}) = ^^^A(r|^^^ (3-2) 
And then we have the equalities 

P(^ - i < L'''^{N,T),R'''^{N,T) <r) (3.3) 
= < X^'^it) <r,l<yi<N,yte [0, T]}) 

/ db/^(-^''-)(T,b|a) 



Now we want to take the limit a, b ^ in (3.2) and a ^ in (3.3) in order to 
prove Theorem 2.1. The following asymptotics in |a| — > for the denominators of 
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(3.2) and (3.3) are given in [15, 17]; as |a| 0, 
/^(r,b|a) 



|b|V(2r) 



(2vr)^/2nf=ir(A;) 
X n {(%-a,)(6,-6,)}x{l + 0(|a|)}, 



l<i<j<N 
N 



yw-A/m N _ _J_ TT r(A:/2) ■v(jv_i)/4 

X n (%-ai)x{l + 0(|a|)}. 

i<i<j<Ar 

By multilinearity of determinants, we have 

/^(-^''■)(r,b|a) 



det 

l<i,j<N 



1 



£ + r 



e + r 



e + r 



N\ \e + r 



N 



T|n|V{2(£+r)2} 



X det 

l<i,j<N 



Sm TT 



det 

l<a,/3<Ar 



sm TT 



where |np = ^^^^ We wiU use the series expansion 



£ + a 

sm I TT- n 

£ + r 



Cp{n){an)^ 



with 



Cp(n) 



(_l)(p-i)/2 



TT 



exp 



n + (-l)P+^exp 



2p! V-^ + ^/l \-^ + ^ / V + 

p e Nq. We will also use the property of the Vandermonde determinant 



det 



n (%-«^)- 

l<i<j<N 
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Then 



/^(-^''•)(T,b|a) 

TV 



(3.9) 



N\ \i + r 



-7r2r|n|2/{2(^+r-)2} 



X 



0<pi<-<pjv 



V det [afl det 

l<i,j<iV L J l<7,(5<Ar 



Pi 



det 

J l<a,P<N 



N 



A/"! V + r y ^ l<7,<5<iV 



X det 

l<a,l3<N 



sm I TT^r- 'T'/? 



e + r 



n ' 

-' l<i<j<N 



ttj — Qi) X 



N 



det 



n 



i<i<j<Ar 
By (3.8), we have 



c^-i{ns)nj ^ 

{l + (^(|a|)} 

J2 e-'^^^-^/{2(^+'-)^>c„_i(n)n"-^ sin [tt 
a,) X {l + e»(|a|)}. 



-n 



nSN 
(_l)(«-2)/2 



-n 



TT 



2{a-l)\ \i + r 

.neN ^ 



a-l 



n H ; n I sm I TT 



2{i + r 



i + r 



i + r 



■n 



^a 1 ( . .2 7ly/^i. ,\ . f i + bp , A 



(_l)(a-2)/2 / ^ 



2(q;-1)! V-^ + ^ 

X ^2 ( 



a-l 



2{e + r) 



£ + r 



n sm TT 



£ + r 



■n 



(_l)(a-3)/2 / ^ 



a-l 



4(q;-1)! V-^ + ^ 

+(-l)° J]]n"-^exp 



rf ^ exp 



TT^T 2 7rv^(2£ + 6^) 



2(£ + r)2 



n 



2(£ + r)2 £ + r 



-n 
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By the identity (2.2), the above is equal to 



_rp-a/2 



2(«+2)/%V2(a _ 1 

Then we have the following asymptotics in |a| ^ 0, 

/^(-^-)(T,b|a) 

_ (_l)^{''V+l)/2 _rj,-N{N+l)/A 



2(£ + r) 



2T 



(3.10) 



X det 

l<a,P<N 



2{i + r) 2i + bf, 



2T 



2T 



+ (-i)"e„-i 



2{i + r) 



2T V2T 



X n (a.-a.)x{l + 0(|a|)}. 



l<i<j<N 



From (3.6) to (3.9), we performed the series expansion (3.7) in only one of the 
two determinants in (3.6). If we perform the series expansions in both of them, we 
have the following estimate. 



/r'^HT,b|a) 



TV! \i + r 



X 



E 



g-7r2r|n|V{2(^+r)2} 



0<P1<---<PN 

- E 

0<(3i<---<(5fjv 
N 



det a' 

l<i,j<N I 



det 

l<a,l3<N 



det 

l<i',j'<N 



det 

J i<Q',^'<Ar 



det 



.neN 



(--\- r J l<a,l3<N 

X n {{aj-a,)ibj-k)]x{l + 0{\a\,\h\)}. 



l<i<j<N 
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By (3.8), we see 



nSN 



a+/3)/2 



TT 



4(a- 1)! \i + r 



a+0-2 



^ g-^2TnV{2(^+r)2}^a+/3-2 



X < exp 



(_l)("+/3)/2 



n) + (-ir + (-l)'' + (-l)"+^exp 

a+/3-2 



-n 



n 



4(a- 1)! V^ + r 



X 



E 



exp 



2(£ + r)2 
+ 5^exp( 



^2 ^ 27r^xAT£ ^ ^ 



2(£ + r)' 



2 ^a+/3-2 



By the identities (2.2), the above is written as 



2(°+/3+i)/2^i/2(c, - l)!(/3 - 1)! 

'2(£ + r) 2£ 



rp-{a+^-l)/2 



X <^ 



^ + {-ife^+0-2 



2{i + r] 



2T 



,0 



^2T V2T 

Therefore, we have the following asymptotics in |a| — >■ and |b| — > 0, 

/^(-^-)(T,b|a) 

_rp-N^/2 



(3.11) 



N 



2^y^7r^/^Yl^^^r{k)^ 



X det 

l<a,/3<7V 



a+0-2 



2{e + r) 2i 



^ ' V2T 



.(-:fe_(^,o 



X 



n {(a,-a,)(6,-6,)}x{l + 0(|a|,|b|)}. 



l<j<j<iV 



By the equalities (3.1) and (3.2), and using the asymptotics formulas (3.4) and 
(3.11), we can take the limit 



lim P f - £ < L^'^iN, T), R^'^iN, T) < r) = lim 

|a|,|b|^0 V V / V / y |^|_n^|_ 



/^(-^■'-^(r,b|a) 

\b\Uo /^(T,b|a) ' 



and then the formula (2.3) is obtained. 
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By the equalities (3.3) and using the asymptotics formulas (3.5) and (3.10), we 
have 

(3.12) 



P{ - e< L°'^(Ar, T), i?°'^(Ar, T) < r) 
/ db/^(-^''-)(T,b|a) 



= lim " ^ 

laKO 



2iv(iv+i)/4 niLir(A;/2) 



/ 



X / (ib det 

+(-i)^e,_i 

(_l)iV(JV+l)/2 

2iv(Jv-i)/4]-[^^^r(A;/2) 



2(£ + r) 2£ 



+ 



/ 



X / dx det 6i_i 



V2T ' ^2T \/2r 
2(£ + r) 

2T ' V^T 



2{i + r) 2i 



2T V2T 



X. 



2T 



, Xj 



If we apply the following lemma, which is known as the de Bruijn identity [5], 
the formula (2.4) is obtained. 

Lemma 3.1 Let Zi{x), 1 < i < N be an integrable piecewise continuous function on 
a region A C M. Let W^(A) = ja; G : G A, 1 < i < ivj. Then 



f 



dx det 

A (A) 1<«J<JV 



Zi{Xj) 



Pf 

l<i,j<N 



(3.13) 



where 



Zi{x)dx, 



I 



det 



(a;i,X2)6A2:i;i<a;2 



Zi[Xi) Zi{X2) 
Zj{xi) Zj{x2) 



dxidx2, 



and 



Zij — 



hj, if l<t,J<N, 

h, ^/ l<t<N,j = N + 1, 

-Ij, if i^N + l,l<j<N, 

0, if i^j^N + 1. 



The proof of Theorem 2.1 was then completed. 
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4 Proof of Theorem 2.2 

For h > 0, consider the restricted region in W^, 



(/i) = jcc = {xi,...,xn) eW%:xN< /i}. 



The transition probabihty density of a one-dimensional Brownian motion in an in- 
terval (0, h), where two absorbing walls are put at x — and x — h > 0, is given 

pTm-) = I E - (T) (T) ■ 

neN 

t>Q,x,ye W%{h). For a, b G W^(/i), let 

/S^'^)(r,b|a)= dct te^)(T,6,|a,)l. 

i<i,j<N L J 

By the Karlin-McGregor formula, it gives the probability density that the absorbing 
Brownian motion in W^(/i) starting from a G W^(/i) at time arrives at b e W^(/i) 
at time T. 

Then we have the equalities, for /i > 

p(if--'-(iV,r)</^) = ^J^^^^^^^^\ a,bewc(/,), (4.1) 

/■ db/^('^)(T,b|a) 
The following asymptotics in |a| — >■ were reported in [13, 15, 19]; as |a| — > 0, 



2iV/2 



b|a) = ^ ^ y-iv(2iv+i)/2^-|b|V(2T) (43) 



TV 



l<i<j<N k=l 



^ ' m=l ^ ' 



N 

X n («i-«?)n«^x{i+^(i^i)}' 

i<i<j<N k=l 

By multilinearity of determinants and by performing the series expansion of a 
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sine function in a determinant, 
/«""(r,b|a) 

N\\h 



X det 



l^sin ^ 



-7r2r|n|2/(2?i2) 





det 






)] 


sin 1 






l<a,l3<N 





_l'^N(N-l)/22N^N^ 



hNiN+1) n-^^ r(2m) 



X det 

l<a,0<N 



X n I 

l<i<j<N 



n sm 



.neN 



2 2 



)n«^^{i+^(i^i)}- 



fe=i 



Here we note that 



T„=/PI.»)„2«-l,i„ 



neN 



V y 



TT 



^ 2 , 7rv^&£ 



-n + 



J16Z 



-r-"e2a-i 



2h hp 



-n 



2T' a/2T ' ' 



where we have used the identity (2.2). Then (4.5) is written as 



2N{N-l)/2^N/2 ji^^^ V{2m) i<«,^<JV 



det 



e 



2a-l 



2h 



2T V2T 



X n I 

l<i<j<7V 



2 2 



AT 

na.x{l + 0(|a|)}. 



k=l 
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On the other hands, in |b| — )■ 0, (4.5) has the asymptotics 



X det 

l<a,l3<N 



\nf-'] det [nf-i] 

L J l<7,5<Af L ^ J 

AT 

X n {{a',-a^){b'^-b^)}l[{a,b,]x{l + 0{\^\,\h\)} 



l<i<j<N 



k=l 



det 



E 



2a+2/3-2 -■ 7r2rn2/(2h2) 



/i^(2^+i) nil=i r(2m)2 i<a;p<N 



l<i<j<Af fc=l 

By using the identities (2.2), the above is written as follows, 
/S^'^(7^,b|a) 

^ I _rp-N{2N+l)/2 



(4.7) 



2N{2N-l)/2^N/2 Yl^^^ r(2 



m] 



X det 

l<a,^<Ar 



e 



2(a+/3-l) 



2h 



,0 



l<i<j<N 



k=l 



Using the asymptotics formulas (4.3) and (4.7), we can take the limit |a|, |b| 
in (4.1). Then (2.5) is obtained. 

By the equality (4.2), and using the asymptotics (4.4) and (4.6), we have 



/ db/5('^)(r,b|a) 



lim 

|a|^0 



rp-N/2 



j 



dh det 

1 



02^-1 



2h bj 



2N(N-l)/2 Yl^^^T{k) 



dx det 

\(h/V2T) ^<M<N 



e 



2T V2T 
2h 



2i-l 



1 -^3 



Then, we use the de Bruijn identity given in Lemma 3.1 and the formula (2.6) 
is obtained. The proof of Theorem 2.2 was then completed. 
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